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ABSTRACT 

Context. In previous works of this series, we have shown that late B- and early A-type stars have genuine bimodal distributions of 
rotational velocities and that late A-type stars lack slow rotators. The distributions of the surface angular velocity ratio fl/Qciii (i\rit is 
the critical angular velocity) have peculiar shapes according to spectral type groups, which can be caused by evolutionary properties. 
Aims. We aim to review the properties of these rotational velocity distributions in some detail as a function of stellar mass and age. 
Methods. We have gathered v sin ; for a sample of 2014 B6- to F2-type stars. We have determined the masses and ages for these objects 
with stellar evolution models. The (reff , log L/Lo)-parameters were determined from the uvby-fi photometry and the HIPPARCOS 
parallaxes. 

Results. The velocity distributions show two regimes that depend on the stellar mass. Stars less massive than 2.5 M© have a unimodal 
equatorial velocity distribution and show a monotonical acceleration with age on the main sequence (MS). Stars more massive have 
a bimodal equatorial velocity distribution. Contrarily to theoretical predictions, the equatorial velocities of stars from about 1.7 Mq 
to 3.2 Mq undergo a strong acceleration in the first third of the MS evolutionary phase, while in the last third of the MS they evolve 
roughly as if there were no angular momentum redistribution in the external stellar layers. The studied stars might start in the ZAMS 
not necessarily as rigid rotators, but with a total angular momentum lower than the critical one of rigid rotators. The stars seem to 
evolve as differential rotators all the way of their MS life span and the variation of the observed rotational velocities proceeds with 
characteristic time scales 6t ^0.2 fms. where tus is the time spent by a star in the MS. 
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1. Introduction 

A great number of phenomena are involved in establishing the 
observed surface rotational velocities of stars. There is a first 
redistribution of angular momentum among the fragments de- 
termined by the fragmentation process of protostellar clouds 
(Bodenheimer, 1981; Matsumoto et al., 1997), which need to be 
dissipated in difl'erent ways before they reach the birth-line. This 
may happen mostly through magnetic braking (Mouschovias 
& Morton, 1985b,a) and bipolar overflows (Pudritz, 1985). 
Following this early evolutionary phase, the amount of angular 
momentum stored by a star before it reaches the zero-age main 
sequence (hereafter ZAMS) can be determined by the gains and 
losses regulated by the accretion, winds (Schatzman, 1962; Shu 
et al., 2000), and locking with the circumstellar environments 
with magnetic fields (Konigl, 1991; Edwards et al., 1993). The 
shaping of the internal distribution of the stored angular momen- 
tum is certainly a matter of hydrodynamical instabilities (Endal 
& Sofia, 1981; Staufi^er et al., 1984; Pinsonneault et al., 1990; 



Send offprint requests to: Frederic Royer, 
e-mail: frederic . royerOobspm . fr 

* Full Tablet is only available in electronic form at CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg . fr/cgi-bin/qcat?J/A+A/ 
** Appendix is only available in electronic form at 
http : / /www . aanda . org 



MacGregor & Brenner, 1991; Keppens et al., 1995; Soderblom 
et al., 1993, 2001), couplings between radiative and convec- 
tive regions (Latour et al., 1981), and by the lockings of layers 
through magnetic fields (Barnes, 2003), whether they be fossil 
or created by various instabilities (Spruit, 1999, 2002). 

If the quoted processes were to lead to establishing only the 
rigid rotation profile in stars on the ZAMS, the distributions of 
angular velocity rates Q/Qcrft (^^crft is the critical rotational ve- 
locity) would probably be unimodal and Maxwellian (Deutsch, 
1970). However, according to differences in the initial chemical 
composition and perhaps other conditions, which may regulate 
the effectiveness of the processes causing the redistribution of 
the angular momentum in the stellar interior, multimodal dis- 
tributions of the Cl/Clciit ratio can be expected. Indeed, bimodal 
rotational velocity distributions are found by Guthrie (1982) for 
late B-type stars in clusters. The bimodality is observed among 
young solar mass stars in Orion (Attridge & Herbst, 1992; Choi 
& Herbst, 1996; Herbst et al., 2001; Barnes, 2003) and is present 
also among the A-F stars of the main sequence (hereafter MS), 
where it is correlated with several spectroscopic characteristics: 
slow rotation favors the appearance of the Am and Ap phenom- 
ena (Abt & Morrell, 1995). However, the slow rotation of Am 
stars can be caused by tidal braking because they are known to 
be binaries (Debemardi, 2000), while the Ap characteristics can 
be favored by the magnetic braking in MS evolutionary phase 
(Hubrig et al., 2000; St?pieh, 2000). 
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Using a highly homogeneous set of v sin ; parameters deter- 
mined for about 11 00 stars, divided into six spectral type groups, 
carefully cleansed of objects presenting the Am/Ap phenomena 
and of all known binaries, Royer et al. (2007) showed that late 
B and early A-type MS stars have genuine bimodal distributions 
of true equatorial rotational velocities. These authors also no- 
ticed a striking lack of slow rotators among the intermediate and 
late A-type stars. They drew attention to the peculiar behavior of 
the frequency distribution of the surface angular velocity ratios 
Q/Qcrit in different spectral type groups. Because these groups 
are characterized by different average ages, the question is raised 
whether this peculiarity could be the consequence of some dif- 
ference in the evolution of the rotation according to the stellar 
mass. 

It is then tempting to review anew in some detail the velocity 
distributions of stars in the MS phase studied above using a bet- 
ter resolution in mass and age than in Royer et al. (2007). This 
may enable us to detect possible signatures in rotational velocity 
properties induced by the stellar formation characteristics and/or 
conditioned by the evolution that follows the ZAMS evolution- 
ary stage. 

The aim of this paper can be summarized as follows: (i) we 
intend to estimate fundamental parameters, mainly masses and 
ages of all stars in the sample, to identify possible statistical in- 
dications on differentiated evolutionary characteristics of the ro- 
tational velocity as a function of the stellar mass and age; (ii) to 
determine what kind of internal rotational distribution may char- 
acterize the stars according to their mass and age during their 
MS phase, which may have some effect on the observed surface 
rotational velocities. 

The selection of the stellar sample studied and the rotational 
velocity data are described in Sect. 2. The determination of 
masses and ages is presented in Sect. 3. The rotational velocity 
distributions are presented in Sect. 4, with details on the statisti- 
cal processing of stellar samples to obtain the equatorial veloc- 
ity distributions from the observed vsin/ values. The evolution 
of these distributions with age is also discussed. Section 5 com- 
pares the observed evolution of rotational velocities with differ- 
ent theoretical models. Finally, the results are summarized and 
discussed in Sect. 6. 



2. Rotational velocity data 

The sample of stars used in the present work is primarily deter- 
mined by the objects studied in Royer et al. (2002a,b, 2007, here- 
after Papers I, II, and III). Because we attempt to study the rota- 
tional characteristics only of A-type stars in the MS evolutionary 
phase, the basis of our sample is the selection of dwarf stars with 
luminosity class from V to IV made in Paper III. However, the 
selection on the spectral class, from B9 to F2, is similar to a cut 
in temperature and strongly affects the distribution in the plane 
mass-age. The first half of the MS, corresponding to young stars, 
is not sampled in the hot end of the selection (spectral types B9 
to A 1). In order to cover the entire age sequence for stars more 
massive than about 2.5 Mq that correspond to early A and late B 
stars, we chose to extend the spectral type range up to B6. 

The final selection of objects that comply with the purpose 
of the present study, however, was made using the masses and 
ages derived using the effective temperatures and bolometric lu- 
minosities determined in the present contribution as detailed in 
Sect. 3. 2. We briefly recall the main characteristics of the basic 
sample gathered in the above cited papers in the following sub- 
sections. 



2.1. vsin i data sources 

Our study is based on four different sources of v sin i data, from 
which we collected both B- and A-type stars: 

- Paper II provides a large and homogeneous v sin / data set 
for A-type stars by merging data from Abt & Morrell (1995) 
with accurate v sin / determined by Fourier transforms (FT). 
The first set is scaled to the later. 

- Abt et al. (2002, hereafter called ALG) provide homoge- 
neous V sin / for 1092 B-type stars derived from the FWHM 
of He 1 447 1 and Mg ii 448 1 A lines calibrated into the v sin / 
scale of Slettebak et al. (1975). This is the extension of the 
work from Abt & Morrell (1995) toward B-stars and it is 
limited to stars brighter than V = 6.5 mag. 

- Levato & Grosso (2004, hereafter called LG) give the south- 
ern counterpart of ALG's work: v sin i for 1027 B-type stars 
determined with the same method. 

- Huang et al. (2010, hereafter called HGM) study rotation in 
a sample of B-stars (in open clusters and field). Their v sin ; 
determination was made by fitting synthetic model profiles 
of He 1 447 1 and Mg n 448 1 A Unes. Only their field star data 
are taken into account in this work. 

This study covers the spectral types from B6 to F2. The num- 
ber of stars per spectral type is given in Table 2. In the total list 
of our sample stars (Table 1) we provide the source of the v sin i 
data: sources 1, 2, and 4 (and their combinations) are taken from 
Paper II, source 8 is taken from ALG, 16 from LG and 32 from 
HGM. 

2.2. Merging and scaling v sin i data 

Although in Paper III a merging with data from ALG was al- 
ready performed to increase the number of stars corresponding 
to the spectral types B9 and B9.5, a new merging was performed 
in this work on a wider range of spectral types. 

The data sets from ALG and LG were first merged, consider- 
ing only stars later than B6. Then the merged sample ALGULG 
was merged with the A-type star sample from Paper II. Finally 

V sin i of field B-type stars from HGM were merged with the full 
sample. The different v sin ; scales are compared in Fig. 1 for the 
94 stars in common between ALG and LG, and for the 199 stars 
in common between Paper II and ALGuLG, and for the 44 stars 
in common between Paper IIUALGuLG and HGM. 

The linear regression lines between the different scales are 
determined with GaussFit (Jefferys et al., 1998a,b), a robust 
least-squares minimization program, to obtain empirical func- 
tions: 

V sin iALG = 1 .28±o.03 v sin j'lg - 8.5±2.5. (1) 

V sin j'paperii = 0.967±o.i9 V sin j'alg + 0.84±o.9. (2) 

V sin j'paperii = 0.971±o.04 V sin j'hgm + 13.7±5.i. (3) 

A 10% error on the v sin i values was assumed for all data when 
using GaussFit to derive these relations. 

Although the v sin ; from Slettebak et al. (1975) are affected 
by the gravitational darkening effect, Howarth (2004) showed 
that they are underestimated compared with more recent grav- 
itational darkening-dependent determinations of v sin i (Fremat 
et al., 2005). Except for the assumptions that underlie the FT 
method, the v sin ; scale defined in Paper II is independent of any 
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Fig. 1. Comparison of v sin i data: (a) Values for the 94 stars (with spectral type equal to, or later than, B6) in common between Abt et al. (2002) 
and Levato & Grosso (2004). The dashed line stands for the one-to-one relation and the solid line is the result from Eq. 1. (b) Values for the 199 
stars (with spectral type equal to, or later than, B6) in common between merged data from Abt et al. (2002) and Levato & Grosso (2004), and 
Royer et al. (2002b). The solid line is the result from Eq. 2. (c) Values for the 44 stars (with spectral type equal to, or later than, B6) in common 
between Huang et al. (2010) and previously merged data. The solid line is the result from Eq. 3. 



other calibration. This scale is moreover confirmed by a more ro- 
bust FT method (Reiners & Royer, 2004) and proves to be con- 
sistent also with the v sin ; determinations by Erspamer & North 
(2003) and Fekel (2003). Recent results from Di'az et al. (201 1) 
provide accurate measurements of v sin i using FT of the cross- 
correlation function for 25 1 A-type stars in the southern hemi- 
sphere. Their method is much less sensitive to blends than the in- 
dividual line measurements made in Paper I, and their v sin i val- 
ues are 5% higher for v sin i > 150 km s ' . This elfect is smaller 
than the 10% error assumed for our merged sample, however, 
and because of the late knowledge of their results, we chose to 
leave our scale unchanged. We therefore used the FT scale of 

V sin i defined in Paper IT and used Eqs. 1, 2, and 3 to scale all 

V sin i data to this scale. The full sample contains 2014 stars and 
the homogenized v sin i are given in Table 1 . 



2.3. Chemically peculiar and binary stars 

In order to have a stellar sample that is as clean as possible 
from objects whose rotation could have been modified by tidal 
or magnetic braking, all known chemically peculiar stars (CP) 
and "close" binary stars (CB) were discarded. The selection was 
made with the same criteria as we used in previous papers. 

CP stars (chemically peculiar stars): the catalog of Ap and 
Am stars from Renson & Manfroid (2009) and the spectral clas- 
sifications given by ALG, LG and HGM were used to identify 
peculiar stars. 

CB stars ("close" binary stars): This category of stars was se- 
lected on criteria based on HIPPARCOS and spectroscopic data. 
Most of the selected stars are in the HIPPARCOS catalog (ESA, 
1997). The binaries detected by the satellite with AHp < 4 mag 
are flagged as CB stars. The Ninth Catalog of Spectroscopic 
Binary Orbits (Pourbaix et al., 2004) was used to complete the 
identification. 

All stars that do not obey CB or CP criteria are simply called 
"normal". The question of identification completeness is studied 
in Paper 11, which concludes that the fraction of CP stars to aU 
stars in the merged sample is fairly constant for aU magnitudes 
brighter than V = 6.5 mag and represents roughly 15% of the 
objects. 



3. Fundamental parameters 

In the present contribution, the rotational properties of stars are 
studied as a function of age and mass following a statistical ap- 
proach. The studied objects were gathered into groups that con- 
tain enough data to warrant that the projection effect in v sin i can 
be statistically corrected. The subtle deviations caused by fast 
rotation efifects, in particular those affecting the v sin ; parameter 
(Zorec et al., 2002; Fremat et al., 2005), can easily be obliter- 
ated in the averaging process by the uncertainties affecting the 
individual fundamental parameters. Therefore we inferred stel- 
lar masses and ages from the apparent effective temperatures and 
bolometric luminosities in the present approach as if the objects 
were at rest. 

3.1. Effective temperatures and bolometric luminosities 

The effective temperatures of stars were derived using uvby-p 
Stromgren photometry and the calibration of Stromgren pho- 
tometric color indices into T^s. The calibrations employed 
are from Moon & Dworetsky (1985), its revisited version for 
some spectral types and luminosity classes obtained by CasteUi 
(1991), and the corrected and extended one for the hottest stars 
by Napiwotzki et al. (1993). The adopted T^ff for a given star 
is the average obtained as follows. The {b - y), mi, c\ and 
P indices were taken from Hauck & Mermilliod (1998) and 
used with their respective imcertainty. The photometric indices 
were used following a Monte Carlo sampling within intervals 
{{X) - IcTx, {X} + 2o-x), where crx is the dispersion of each in- 
dividual photometric index X. For each star we obtained 6500 
reff-determinations, whose average and the corresponding dis- 
persion are the adopted T^g ± crj^ values presented in Table 1. 
In Fig. 2 we show the obtained effective temperatures against 
the photometric colors {b - y) corrected for interstellar extinc- 
tion. For some stars in the sample of Huang et al. (2010) uvby-ji 
photometry is not available, so that we used the T^.^ estimates 
derived by these authors. 

Sometimes, more or less systematic deviations can be ex- 
pected between the effective temperatures determined from pho- 
tometric indices with those based on bolometric fluxes (Zorec 
et al., 2009). While the details of the stellar energy distribution 
can be important when T^s is determined with photometric in- 
dices, they have a marginal incidence in the refif-determination 
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Table 1. (extract) List of the 2014 B6- to F2-type stars, with their spectral type and fundamental parameters. 



HD 


Spec, type 


log Teff 




logL/Lo 






0"M 


t/tus 


""'/'MS 


vsm( 

(kms-i) 


Source 


Classification 


3 


AlVn 


3.957 


0.008 


1.6387 


0.0725 


2.36 


0.08 


0.685 


0.075 


228 


4 


CP 


203 


F2IV 


3.832 


0.003 


0.6669 


0.0146 


1.42 


0.01 


0.449 


0.048 


170 


4 




256 


A2IVA' 


3.943 


0.009 


2.0664 


0.1160 


2.79 


0.12 


0.983 


0.040 


241 


5 




319 


AlV 


3.945 


0.009 


1.4761 


0.0306 


2.18 


0.03 


0.592 


0.062 


59 


5 


CP 


431 


AVIV 


3.889 


0.011 


1.5825 


0.0554 


2.23 


0.06 


0.899 


0.042 


97 


4 


CB 


560 


B9V 


4.035 


0.004 


1.8414 


0.0293 


2.74 


0.03 


0.417 


0.054 


249 


1 




565 


A6V 


3.906 


0.004 


1.6779 


0.0422 


2.35 


0.05 


0.900 


0.025 


149 


1 




584 


B7IV 


4.116 


0.003 


2.1390 


0.0740 


3.38 


0.10 


0.140 


0.127 


12 


8 




709 


B8 


4.078 


0.002 


3.2140 


0.5320 


5.17 


1.60 


0.988 


0.060 


264 


32 


CB 



Notes. (Tj and cr^ are the errors in logarithmic T^-g and logarithmic luminosity. The homogenized v sin ( is given (with corresponding source) as 
well as the classification (CP, CB, blank stands for normal). The source flag is the same as in Paper II, and additionally 8 for ALG, 16 for LG, and 
32 for HGM. 
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Fig. 2. Color-effective temperature diagram for the single normal stars 
selected for the present study. 
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Fig. 3. Comparison of effective temperatures calculated in this paper 
from the uvby - p photometry and those for 135 stars obtained using 
bolometric fluxes. The dashed line is the one-to-one relation. 



with bolometric fluxes (Fremat & Zorec, 2003; Fremat et al., 
2003). We compared the effective temperatures of 135 stars 
obtained in this work with those determined with bolometric 
fluxes. Techniques for determining eff'ective temperatures with 
integrated bolometric fluxes are given in Zorec et al. (2009) and 
they are presently used to determine the effective temperature of 



600 A- andF-type stars (in preparation). The comparison (Fig. 3) 

shows a great consistency. 

The bolometric luminosity of stars is estimated with 



Mboi = Mv + BC(T^g), 



(4) 



where My is the visual absolute magnitude and BC(Teff) is the 
bolometric correction (Lang, 1992). To calculate My we used 
(i) the apparent visual magnitude of the UBV photometry and 
the corresponding uncertainty; (ii) the HIPPARCOS parallaxes 
recently re-determined by van Leeuwen (2007) and taking into 
account the respective uncertainties; (iii) the interstellar color 
excess E(B - V) estimated using the uvby-fi photometry and 
the calibration of intrinsic colors by Moon & Dworetsky (1985). 
Each adopted Mboi magnitude is the average of all determina- 
tions by Eq. 4 following a Monte Carlo sampling of the un- 
certainties affecting T^s and the parameters entering Pogson's 
relation for My. The adopted bolometric luminosity parameter 
log L/Lq, given in Table 1 with its Icr uncertainty, was derived 
by adopting = 4.742 mag. 



3.2. Stellar masses and ages 

To infer the individual stellar masses and ages, we used the pairs 
(Feff, log L/L0) obtained above with their respective 2cr devia- 
tions as entries in the evolutionary diagrams of non-rotating stars 
calculated by Schaller et al. (1992) that we display in Fig. 4. The 
interpolation of masses and ages also proceeds with a Monte 
Carlo sampling of the (reft, log L/L©) entries. Table 1 gives the 
average (M/Mq, f/fMs), i-e. mass and fractional age parameters 
with the Icr dispersion that are the results from each interpola- 
tion procedure. We use here the notation t for the age in years, 
and t/tus is the life span from the ZAMS to the TAMS (termi- 
nal age main sequence) of a star with the given mass M/Mq. 
Figure 5 shows the distribution in the mass-age diagram of the 
MS stars used in the present study. The vertical cuts in effective 
temperature in Fig. 4 that appear because of the selection based 
on MK spectral type are reflected in Fig. 5 by the limits in mass 
varying with age. Table 2 gives some statistical estimators that 
describe the characteristics of the stellar mass distribution tabu- 
lated against the spectral types of the present analysis. 

The stars with large error bars on mass and/or age (ct/Mq > 
0.3, (Tt/tff^ > 0.15) were discarded for the next steps of this work. 
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Table 2. Parameters of the mass distributions in different spectral type 
intervals. 



Fig. 4. HR diagram of the studied stars. Only the single normal stars 
in the shaded region were selected for the present study. Filled circles 
are for stars in the sample studied by Royer et al. (2007). Open circles, 
open diamonds and asterisks stand for complementary stars taken from 
Abt et al. (2002), Levato & Grosso (2004) and Huang et al. (2010) re- 
spectively. The evolutionary tracks are from Schaller et al. (1992). 




Fig. 5. Distribution in the {t/tus, M/Mq) plane of the data concerning 
the single normal stars selected for the present study. The error bars 
in both coordinates are indicated in gray. The symbols have the same 
meaning as in Fig. 4. 



Spectral 


# Ml cr. 


M2 (To 


M3 


type 


stars (Mo) (Mq) 


(Mo) (Mo) 


(Mo) 



B6- 


-B7 


135 


4.02 


0.51 


4.24 


0.64 


4.09 


B8 




182 


3.43 


0.54 


3.52 


0.54 


3.44 


B9- 


-B9.5 


328 


2.99 


0.45 


3.08 


0.55 


2.99 


AO 




193 


2.57 


0.36 


2.66 


0.42 


2.59 


Al 




193 


2.51 


0.38 


2.59 


0.35 


2.53 


A2 




194 


2.47 


0.43 


2.50 


0.41 


2.50 


A3 




170 


2.27 


0.38 


2.39 


0.45 


2.30 


A4- 


-A6 


179 


2.11 


0.31 


2.24 


0.43 


2.14 


A7- 


-A9 


180 


1.93 


0.31 


1.99 


0.49 


1.93 


FO- 


F2 


194 


1.68 


0.26 


1.84 


0.39 


1.74 



Notes. The values Mi and cri are the center and sigma of the Gaussian 
fit of the distribution, M2 and cr2 are the mean and the standard deviation 
of the distribution and M3 its median value. 

Table 3. Mass and age estimates using evolutionary tracks with rotation 
for a test star with apparent rotation-less mass M = 3 Mo. 





M(fi)/Mo 




m)/tusm 






t/tus = 0.2 0.5 


0.8 


///MS = 0.2 


0.5 


0.8 


0.5 


3.05 2.91 


2.83 


0.20 


0.49 


0.78 


0.8 


3.14 3.00 


2.85 


0.19 


0.46 


0.70 


0.9 


3.18 3.07 


2.89 


0.19 


0.44 


0.67 


1.0 


3.22 3.13 


2.93 


0.18 


0.42 


0.63 



3.2.1. Uncertainties regarding tine stellar masses and ages 

Because models of stellar evolution with rotation can produce 
evolutionary tracks that can differ significantly from those with- 
out rotation (Maeder & Meynet, 2000; Ekstrom et al., 2008, 
2011), we may wonder whether the masses and ages obtained 
for our sample stars, which can be intrinsic rapid rotators, are 
still realistic. Because we do not know the actual internal rota- 
tion of a given star in advance, whatever uncertainty we can es- 
timate will necessarily be model-dependent. We can then recall 
that the most frequently used models (c.f. Maeder & Meynet, 
2000; Ekstrom et al., 2008, 2011) are based on two strong as- 
sumptions: (i) all stars initiate their evolution in the MS as rigid 
rotators, which means the total rotational energy stored in the 
object is strongly limited (see discussion in Sect. 5.4); (ii) the re- 
distribution of angular momentum proceeds only over barotropic 
surfaces and produces only "shellular" internal rotation laws. We 
have calculated the masses and rotation-dependent masses and 
fractional ages t{Q.)/tMsi^) suggested by these models for a test 
star that according to rotation-less models has a mass M - 3 Mq 
and lies at evolutionary phases f//MS - 0-2, 0.5 and 0.8. The 
obtained results are given in Table 3. From this table we can 
conclude that the parameters based on rotation-less models used 
in the present work (M/Mq, f(Q)/f/fMs) may slightly underesti- 
mate stellar masses in the first half of the MS, and that the differ- 
ences are larger for a higher ratio Q./Q.cnt- Masses can be slightly 
overestimated for stars in the second half of the MS. Because the 
actual ratio Q/Qcdt is not known, the stellar mass and its uncer- 
tainty are on average (M) - 3.01 ±0.12 Mq. In Table 3 we note 
that rotation-dependent ratios can be systematically lower than 
those derived from rotation-less models and the differences are 
larger the closer the stars are to the TAMS and the higher the ra- 
tio The effect of neglecting rotation on the derivation of 
M/Mq and f/fMS for other stellar masses can be seen in Fig. 6 of 
Zorec et al. (2005), where the authors studied Be stars, which are 
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considered as paradigms for rapidly rotating MS stars. We note, 
however, that the fraction of program stars in the present work 
with Q/Qcrit ~ 1-0 must probably be very low, because the frac- 
tion of Be stars among the B-type star population is low (Zorec 
& Briot, 1997). On the other hand, as soon as Q/Qcrit SO.9, age 
differences become 6it(Q)/tMs) ^0.1, which are on the same 
order or smaller than the age-bins used to estimate the averages 
(v) of true equatorial velocities in the following sections. 

3.2.2. Determining fundamental parameters with 
rotation-dependent evolutionary tracks 

The use of evolutionary tracks with rotation implies that we have 
an estimate of the true equatorial velocity v of the studied star. 
These tracks are given in terms of bolometric luminosities and 
effective temperatures averaged over the rotationally deformed 
object as a function of age and rotational velocities in the ZAMS. 
Accordingly, the observed apparent aspect angle-dependent ef- 
fective temperature and bolometric luminosity must be corrected 
for rotational effects to derive their counterparts for the corre- 
sponding rotation-less object (Fremat et al., 2005; Zorec et al., 
2005). In turn, these must be averaged over a rotationally de- 
formed star, whose degree of geometrical deformation is mea- 
sured by the ratio v/vcrit(M, t), where VctitiM, t) is the stellar crit- 
ical equatorial velocity. Using the curves representing the evo- 
lution of the ratio v(M, t)/vcnt(M, t), the stellar velocity ratio in 
the ZAMS must be identified. Only evolutionary tracks for this 
specific initial velocity ratio, and for several masses around the 
sought one, can finally be used to interpolate the mass and age 
of the studied object. This operation implies a long iteration pro- 
cedure, where the stellar mass M/Mq, its age f/fMS and the incli- 
nation angle / of the rotational axis are iterated simultaneously. 
Knowing that on the one hand only a small fraction of program 
stars may be very rapid rotators, and on the other hand this iter- 
ation will not bring significant improvements in the estimate of 
masses and ages of the remaining stars, the approach used in the 
present work seems justified. 

3.3. The V sin ; parameters and the rapid rotation 

The estimates of v sin i parameters can also be affected by the 
rapid rotation. Rapid rotation induces temperature and gravity 
inhomogeneities in the stellar surface known as the "gravity 
darkening effect". Owing to the lowering induced on the local 
effective temperatures in the equator, these regions do not con- 
tribute efficiently to the Doppler rotational broadening of spec- 
tral lines. The apparent v sin i of rapid rotators must then be cor- 
rected for this underestimation (Fremat et al., 2005). This correc- 
tion can increase our estimates of v/vcnt only if it can be demon- 
strated that in each bin used to average the apparent v sin i param- 
eters, the number of rapid rotators with actual ratios Q/Qc ^ 0.8 
represents a significant fraction of objects. 



4. Distributions of rotational velocities 

In Paper III, distributions of rotational velocities were obtained 
by dividing the stellar sample into groups of spectral types and 
mixing stars of all ages in the MS life span. However, the spec- 
tral type of a star changes during the MS hfe span, so that in a 
given group objects with different masses and ages can display 
roughly the same apparent spectral type. To have a better repre- 
sentation of the percentage of stars as a function of the rotational 
velocity per mass-groups and study the distribution of rotational 



velocities as a function of the stellar mass, the sample of normal 
stars defined above was divided into overlaping mass intervals 
that run from 1.6 to 3.85 Mq. The limits of the intervals together 
with the corresponding number of stars are given in Table 4. For 
each subset, the distribution of vsin; was processed in a way 
similar as in Paper III: 

- the histogram was smoothed using a Gaussian kernel es- 
timator and smoothing parameter h defined by Bowman 
& Azzahni (1997)', and variabihty bands were estimated 
(Bowman & Azzalini, 1997), 

- the smoothed distribution was rectified from the projection 
effect, assuming randomly oriented rotation axes, with the 
Lucy-Richardson method (Lucy, 1974; Richardson, 1972), 

- the rectified distributions were fitted by a sum of two 
Maxwellians to derive the position of the mode(s) and the 
proportion of fast and slow rotators if a significant bimodal- 
ity is present. The Maxwellian distribution used to fit the fast 
part of the velocity distribution has an additional lag param- 
eter {, to be able to shift it toward higher velocities^. 

The distributions are shown for some of the mass intervals in 
Fig. 6, and the results of the fits are given in Table 4. An overview 
of the distribution of the rotational velocities as a function of 
mass for the stellar sample in the present work is displayed in 
Fig. 7. In this figure the full MS life span is considered: < 
t/tus ^ 1- The color scale represents the density in the one- 
dimensional normalized distributions: the bluer the region, the 
smaller the number of stars and vice versa toward the red scale. 

We clearly notice in Fig. 7 two main groups of stars as a func- 
tion of stellar mass with a transition occurring at M » 2.5 Mq. 
In the low-mass part of the diagram, the distribution is unimodal 
and there is a striking lack of objects in the $ v $ 100 km s"' 
interval of velocities. This gap in the distribution in the plane 
(v, M) is also noticeable, although less pronounced, in (v sin i, 
spectral type) as shown in Royer (2009). With the aim of testing 
the effect of our selection of normal stars on the resulting dis- 
tributions, the same process was applied to all stars (normal, CP 
and CB). Table 5 gives the results of the fits for the same mass 
intervals. The lack is still present even when CP and CB stars are 
not removed from the sample, although some minor overdensi- 
ties occur for v 40 km s ' . The effect of discarding CP and CB 
stars cannot account for the lack of slow rotators, and we con- 
clude that stars with low rotational velocities were not created in 
the space volume scrutinized with our data. 

For stars with M ^ 2.5 Mq, the bimodality of the velocity 
distribution distinctly appears. The proportion of the slow rota- 
tors increases to 20% around M » 2.8 Mq and decreases toward 
higher masses. The selection of normal stars removes a large 
proportion of slow rotators (0 ^ v ^ 100 km s"') in this mass 
range, as seen when comparing the percentage of the slow ro- 
tator Maxwellian distributions from Tables 4 and 5, but the bi- 
modaUty remains significant. 

The fast rotator mode of the equatorial velocity distribution 
for normal stars, fi,, behaves monotonically over the full range 



' This smoothing parameter is larger than the estimation using the 
definition by Sheather & Jones (1991), made in Paper III. This is the rea- 
son why the distributions are smoother (Silverman, 1986). This smooth- 
ing parameter has been used in the different parts of the study, both in 
one and two dimensions (see Sect. 5.2). 

^ The lag parameter { allows us to shift the Maxwellian distribution 
along the x-axis: f{x) = (x- if ^fljnlv' exp (- {x - {f / (Sa^)). 
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Fig. 6. Distributions of rotational velocities for normal stars in different mass subsamples: shaded histograms are the observed vsini; solid thick 
lines are the distributions of equatorial velocities and the gray strips are their associated variability bands. The dashed lines stand for the Maxwellian 
fit, whose parameters are given in Table 4. The mass range for each subsample is indicated in the corresponding panel, together with the number 
of stars. The histograms are normalized to fit the probability density scale. 



probability 
density 




V (km s ') 



Fig. 7. Distribution of true rotational velocities v for normal stars as a function of the stellar mass M/Mq. The color scale represents the density in 
the one-dimensional normalized distributions. The bluer the region, the smaller the number of stars and vice versa for the red colored scale. White 
solid lines follow iso-density contours. 
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Table 4. Parameters of the ID velocity distributions for normal stars as a function of the stellar mass. 



Mass 


# 






Slow rotators 








Fast rotators 






range 


stars 


h 


% 




K 


dispersion 


% 




e 




K 


dispersion 


(Mo) 








(kms-') 


(kms-') 


(kms-') 




(kms-') 


(kms-') 


(kms-') 


(kms-') 


(kms-') 


* 1.60-2.00 


191 


0.200 


— 








— 


100 


114±1 


46 ±1 


160 


148 


54 


1.70-2.10 


216 


0.194 


— 








— 


100 


124+1 


46 ±1 


170 


151 


59 


1.80-2.20 


220 


0.177 


— 








— 


100 


128 ±1 


48+1 


176 


160 


61 


* 1.90-2.30 


249 


0.190 


— 








— 


100 


133 ±1 


44 ±1 


177 


169 


63 


2.00-2.40 


276 


0.201 


— 








— 


100 


140 ±1 


45 ±1 


185 


181 


67 


2.10-2.50 


274 


0.227 


3 


25 ±1 


22 


12 


97 


154+1 


27 ±1 


181 


181 


73 


*2.20-2.60 


273 


0.247 


8 


37 ±1 


43 


18 


92 


164 ±1 


18 ±1 


183 


184 


78 


2.30-2.70 


248 


0.260 


12 


44 ±1 


43 


21 


88 


179 ±1 


12 ±2 


191 


187 


85 


2.35-2.85 


282 


0.259 


12 


47 ±1 


49 


22 


88 


210±1 


0±1 


210 


178 


100 


*2.45~2.95 


274 


0.254 


18 


52 ±1 


46 


25 


82 


185 ±3 


14 ±4 


199 


193 


88 


2.55-3.05 


235 


0.254 


20 


61 ±1 


52 


29 


80 


191 ±2 


20 ±3 


211 


214 


91 


2.65-3.15 


209 


0.253 


15 


55 ±1 


46 


26 


85 


183 ±2 


24 ±3 


208 


205 


87 


*2.65-3.35 


259 


0.249 


14 


52+1 


43 


25 


86 


187 ±2 


29 ±3 


216 


214 


89 


2.75-3.45 


225 


0.270 


13 


47 ±1 


40 


22 


87 


201 ±2 


18 ±3 


219 


217 


96 


2.85-3.55 


190 


0.278 


10 


39 ±1 


37 


19 


90 


202 ±2 


20 ±2 


222 


217 


96 


*2.95-3.65 


153 


0.290 


6 


32 ±1 


31 


15 


94 


224 ±2 


Oil 


224 


220 


107 


3.05-3.75 


131 


0.301 


5 


31 ±1 


31 


15 


95 


229 ±1 


0±1 


229 


190 


109 


3.15-3.85 


122 


0.314 


8 


32 ±1 


28 


15 


92 


224 ±2 


15 ±2 


238 


220 


107 



Notes. The defined mass ranges and the corresponding number of normal stars are given. The smoothing parameter h is given for each subsample, 
it is expressed in logarithmic v sin i. For each distribution, the parameters of the Maxwellian fit are listed: percentage of slow and fast velocity 
distributions, mode and dispersion. For slow rotators, the distribution mode jUg is derived from the Maxwelhan fit. The fast rotator distribution is a 
lagged Maxwellian (see text and footnote 2), and its mode is yUr = y/2a+{, where a is the parameter of the Maxwellian distribution. The estimators 
fi' are derived as the position of the maximum directly on the ID distributions. The distributions are shown in Fig. 6 for six of the eighteen mass 
ranges, indicated by asterisks. 



Table 5. Parameters of the ID velocity distributions for all stars (normal, CP and CB) as a function of the stellar mass. The columns are the same 
as in Table 4. 



Mass 


# 






Slow rotators 










Fast rotators 






range 


stars 


h 


% 


Ms 




dispersion 


% 


V2a 


e 




K 


dispersion 










(kms-i) 


(kms-i) 


(kms-i) 




(kms-i) 


(kms" 


) (kms-i) 


(kms-i) 


(kms-i) 


1.60-2.00 


283 


0.196 










100 


122±1 


40±1 


162 


145 


58 


1.70-2.10 


321 


0.189 










100 


127±1 


41±1 


168 


151 


60 


1.80-2.20 


342 


0.174 










100 


132±1 


41±1 


172 


163 


63 


1.90-2.30 


385 


0.194 










100 


142±1 


31±1 


173 


175 


67 


2.00-2.40 


429 


0.207 


5 


50 ±2 


58 


24 


95 


149±1 


32±1 


181 


184 


71 


2.10-2.50 


418 


0.225 


4 


34 ±1 


58 


16 


96 


168±1 


9±1 


177 


184 


80 


2.20-2.60 


426 


0.244 


10 


39 ±1 


37 


19 


90 


176±1 


4±1 


180 


184 


84 


2.30-2.70 


391 


0.263 


15 


40 ±1 


34 


19 


85 


201±1 


Oil 


201 


187 


96 


2.35-2.85 


455 


0.261 


19 


43 ±1 


34 


20 


81 


189±1 


Oil 


189 


181 


90 


2.45-2.95 


447 


0.263 


22 


44 ±1 


37 


21 


78 


192±1 


Oil 


192 


190 


92 


2.55-3.05 


385 


0.275 


23 


45 ±1 


34 


22 


77 


196±1 


Oil 


196 


202 


93 


2.65-3.15 


340 


0.275 


20 


44 ±1 


34 


21 


80 


190±1 


Oil 


190 


193 


90 


2.65-3.35 


415 


0.274 


20 


42 ±1 


31 


20 


80 


203±1 


Oil 


203 


205 


97 


2.75-3.45 


352 


0.287 


19 


41 ±1 


31 


20 


81 


203±1 


Oil 


203 


208 


97 


2.85-3.55 


299 


0.301 


16 


38 ±1 


31 


18 


84 


212±1 


Oil 


212 


205 


101 


2.95-3.65 


233 


0.316 


16 


40 ±1 


31 


19 


84 


211±1 


Oil 


211 


205 


100 


3.05-3.75 


209 


0.337 


16 


36 ±1 


31 


17 


84 


224±1 


Oil 


224 


178 


107 


3.15-3.85 


189 


0.359 


18 


35 ±1 


28 


17 


82 


229±1 


Oil 


229 


208 


109 



of mass and smoothly increases with stellar mass. This variation 
can be fitted by a linear relation in the mass range 1.6-3.5 M©: 

Hr « 41±i.9 M/Mq + 90±5 kms"^ (5) 

4.1. Normal versus peculiar stellar populations 

It was argued by Abt (2009) that slowly rotating A0-A3 "nor- 
mal" stars are stars that did not have enough time to become Ap 
or Am stars yet, and that Ap(SrCrEu) stars need about half of 
their MS lifetime to show their chemical peculiarity. This effect 



would bias the distributions of rotational velocities when one 
considers only "normal" stars and applies a selection by age. 

The averaged age on the MS as a function of v sin i was de- 
rived for different mass intervals in our sample, considering nor- 
mal stars and CP stars separately (Fig. 8). In the first two mass 
intervals (2 < M/Mq < 2.6) there is a significant and system- 
atic shift between normal and CP stars for v sin i < 100 km s~' , 
CP being on average older by 10% of ?ms- One should bear in 
mind that when deriving the fundamental parameters of our tar- 
gets, normal and CP stars are treated in the same way. Therefore, 
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Fig. 8. Mean age (t/tMs) per bin of vsin; for different mass intervals 
(taken from Table 4): (a) 2 < M/Mq < 2.4; (b) 2.2 < M/Mq < 2.6; 
(c) 2.35 < M/Mo < 2.85; (d) 2.55 < M/Mq < 3.05. Filled circles 
stand for normal stars, whereas open squares represent CP stars. The 
bins in v sin / are 40 km s"' wide. Solid lines are the results of linear fits 
on normal stars with v sin / < 110 km s"' and dashed lines the results of 
the fits for CP stars with v sin i > 90 km s ' . 



the bolometric luminosity of CP stars can be overestimated. 
Landstreet et al. (2007) indeed show that the bolometric correc- 
tion for Ap stars is smaller than the one caUbrated by Lanz & 
Catala (1992) by an amount of about 0.2 mag. Their age derived 
from their position in the HR diagram may consequently be also 
overestimated. 

On the other hand, a trend is seen in the variation of the mean 
age of normal stars in the first mass interval (2 < M/Mq < 2.4). 
It increases monotonically from 53% to 68% of tus between 
vsinj = 50 and 110kms"\ and remains constant for faster ro- 
tators. This trend is less pronounced in the next mass bin, and 
completely disappears for more massive stars. This trend agrees 
with that postulated by Abt (2009). Moreover, the variation of 
the mean age with rotation is recovered using v sin i, therefore 
the projection effect can lessen the slope. But the selection here 
was made using mass and age. According to the evolutionary 
tracks in Fig. 4, the A0-A3 spectral types will correspond to dif- 
ferent stellar masses if stars are close to the ZAMS or near the 
TAMS. 

This trend does not seem to explain the bimodahty observed 
in Figs 6 and 7. According to Table 4, the bimodality reaches its 
maximum (in terms of percentage of the full distribution) in the 
mass bin 2.55 < M/Mq < 3.05, for which no obvious trend is 
seen in Fig. 8.d. 



4.2. Evolution of the distributions of rotational velocities 

In the distributions shown in Fig. 6 stars are not distinguished by 
their ages. Below, stellar ages are taken into account in two dif- 
ferent ways: (i) by studying the changes of the rotation velocity 
distributions as a function of time (this section); (ii) by study- 
ing the evolution from the ZAMS to the TAMS of the average 
rotational velocities per mass (Sect. 5). 

The evolution of the distributions of rotational velocity dur- 
ing the MS life time is studied separately for the two groups 
of stars identified in Fig. 7. The cut in mass was chosen to be 
M > 2.4Mq, however, to gather enough stars in the high-mass 
part. We obtained distributions of v for each mass group by sep- 
arating the stars into three age intervals (Fig. 9). They were cho- 
sen so that the number of stars in each of them enabled us to 
rectify the distributions from the projection angle effect. As seen 
in Fig. 5, the young part of the diagram is poorly sampled. Hence 
the first bin is large, < t/tus < 0.5, to warrant statistical reli- 
ability. The two mass bins behave quite differently as a function 
of age. Table 6 gives the parameters that characterize the sub- 
samples and their velocity distributions. 

The low-mass star v distribution does not vary significantly 
in shape with age, and remains unimodal from ZAMS to TAMS. 
It shows an acceleration, however. The mode of the distribu- 
tion increases from about 150km s"' for < t/tMS < 0.5, to 
170 km s ' for 0.5 < t/tus < 0.8 and to about 185 km s"^ for 
0.8 < t/tus < 1. 

The high mass v distribution does exhibit significant changes 
as a function of age. The bimodality is very clear among stars 
younger than t/tus < 0.8, whereas the two peaks blend near the 
TAMS (0.8 < t/tus < !)• Slowly rotating stars are present at any 
age, but the fast rotating part of the distribution decelerates with 
increasing age. In terms of positions of the mode, the distribu- 
tion peaks at about 240 km s"' for the younger half (t/tus < 0.5), 
and stagnates at about 190kms"' for the older part. In terms 
of median of the full distribution, however, the deceleration 
is continuous and the median v decreases from 250kms~' for 
< t/tus < 0.5, to 210 km s"^ for 0.5 < t/tus < 0.8 and to 
185 km s"' for 0.8 < t/tus < 1- 

Using these distributions, the ratio of the average velocity 
of the fast rotating component (derived from the formula of the 
Maxwellian distribution) over the average velocity of the full 
component was derived. This ratio is characterized by the fac- 
tor Kr, which allowed us to correct in Sect. 5.2 the mean values 
of rotational velocities in mixed samples encompassing slowly 
and rapidly rotating stars. It corrects the curves of rotational ve- 
locities against stellar ages for the effect carried by the slowly 
rotating stars, so as to derive distributions and related modes for 
the rapidly rotating stars only. 



5. Evolution of the rotational velocities compared 
with models 

In a first attempt, we tried to apprehend possible dependencies 
of the evolution with age of rotational velocities for different 
masses through the v/Vcrit velocity ratio, where Vcrit is the time- 
and mass-dependent critical equatorial rotational velocity. In a 
second approach, we studied the v/v'zams velocity ratio, where 
vzAMS is the equatorial rotational velocity on the ZAMS. 

The velocity ratios v/vcrit and v/vzams were compared with 
model predictions for three different regimes of internal redistri- 
bution of the angular momentum. The first two predictions are 
for extreme redistributions and are detailed in the present dis- 
cussion. The third type of models were calculated by Ekstrom 
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Fig. 9. Distributions of rotational velocities as a 
function of age for two mass subsamples. Left 
panels correspond to the low-mass bin: 1.6- 
2.4 Mq and right panels correspond to the high- 
mass bin: 2.4-3.85 Mq. As in Fig. 6, shaded 
histograms are the observed v sin solid thick 
lines are the distributions of equatorial veloc- 
ities V and the gray strips are their associated 
variability bands. The mass range for each sub- 
sample is indicated in the corresponding panel 
together with the number of stars, and below, 
the age interval is indicated and goes from 
young stars (0 < f/fms < 0.5) in upper panels 
to stars close to the TAMS (0.8 < f/fMS < 1) in 
lower panels. The histograms are normalized to 
fit the probability density scale. 



Table 6. Parameters of the ID velocity distributions determined in the mass bins and age bins shown in Fig 9. 



Mass 


f/^MS 


# 




Slow rotators 






Fast rotators 








range 


range 


stars 


h 




dispersion 


% 






K 


dispersion 




(Mo) 








(kms-') 


(kms-') 




(km s-') 


(kms"') (kms"') 


(kms-') 


(kms-') 




1.60-2.40 


0.0-0.5 


134 


0.228 






100 


138±1 


26±1 164 


148 


65 


1.005 


1.60-2.40 


0.5-0.8 


217 


0.204 






100 


127±1 


52±1 179 


172 


61 


0.979 


1.60-2.40 


0.8-1 


108 


0.236 






100 


141±1 


40±1 181 


184 


67 


0.995 


2.40-3.85 


0.0-0.5 


84 


0.317 


47 ±1 


22 


89 


221±3 


24±3 245 


241 


105 


1.101 


2.40-3.85 


0.5-0.8 


188 


0.268 


49 ±1 


23 


88 


193±2 


15±3 208 


196 


92 


1.103 


2.40-3.85 


0.8-1 


215 


0.276 


43 ±1 


20 


88 


184±1 


0±1 184 


184 


88 


1.053 



Notes. As in Table 4, the corresponding number of normal stars and the smoothing parameter h (expressed in logarithmic v sin ;) are given for each 
subsample. For each distribution, the parameters of the Maxwellian fit are listed: percentage of fast velocity distributions, mode and dispersion. 
The mode of the slow rotator distribution /i^ is derived from the Maxwellian fit. The fast rotator distribution is a lagged Maxwellian (see text and 
footnote 2) and its mode is /ij = y/la + C. The estimator fi[. was derived as the position of the maximum directly on the ID fast distributions. 
Correction factors were derived as = (2 a -s/l/n + £)/{v), where (v) is the mean of the full distribution. 




100 200 300 400 500 
rotational velocity (km s"') 



100 200 300 400 500 
rotational velocity (km s"') 



et al. (2008) for stars with 3 Mq, the only common mass with the 
present paper 

The ZAMS was considered as the initial state of the MS 
stellar evolutionary phase. Nevertheless, as a consequence of all 
possible physical conditions prevailing during the evolutionary 
stages preceding the ZAMS, a given star on the ZAMS could 
be either a rigid rotator or a differential rotator. If the star is a 
differential rotator, the total angular momentum stored in the ob- 
ject can be higher or lower than the limiting value supported by 
a critical rigid rotator. The characteristics of the stellar struc- 
ture depend on the amount of rotational energy stored and on 
the internal distribution of the angular velocity. The competition 
between gravitational and centrifugal forces determines the final 
stellar internal structure, the meridional circulation, and the mix- 
ing phenomena of chemical elements that put constraints to the 
subsequent stellar evolution in the MS (cf. Maeder & Meynet, 
2000). 

The mass-loss rates inferred for A-type stars are on the or- 
der of 10"'^ to 10"'"Mayr-i (Lanz & Catala, 1992; Lamers 
& Cassinelli, 1999), which for evolutionary time scales rang- 
ing from 10^ to some 10"^ yr implies that the amount of angular 



momentum removed in these objects during their MS life span is 
probably negligible on average. We assumed then that all studied 
stars evolve during the MS phase conserving their total angular 
momentum. 

The amount of rotational kinetic energy stored by a star is 
currently estimated using the ratio t = A'/IWI, where K is the 
kinetic energy and W is the gravitational potential energy. At 
rotational energies approaching t » 0.27, the object is dynam- 
ically unstable. For r x 0.14 it becomes secularly unstable, in 
the sense that its axial symmetry can be broken and transformed 
into a three-axially symmetric body (Tassoul, 1978). If at any 
time during the MS evolution a star has a total angular momen- 
tum higher than permitted by the critical rigid rotation, the star is 
said to evolve as a neat differential rotator. We refer to the evolu- 
tion at low rotational energy regime if stars have r < t"^^ , oth- 
erwise they will be said to evolve at high rotational energy. The 
notation t"^^ stands for the critical, or upper possible energy 

born by a rigid rotator, which ranges from t < t"^^"' - 0.0041 
to 0.0053 for stars with masses from M = 1.5 to 3Mq, respec- 
tively. Assuming that stars are permanently axisymmetric, the 
evolution at high-energy regime implies that <t<0.14. 
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In the present discussion we assumed that stars evolve in the 
MS at low rotational energy regime, i.e. they have r < t^^^. 
Moreover, in our calculations we still limited the total angular 
momentum to / < Jj^^^, where for a given mass, Jj^^^ cor- 
responds to the angular momentum supported by a critical rigid 
rotator on the TAMS. This choice is mandatory for stars that are 
assumed to be rigid rotators on the ZAMS and they evolve in the 
MS by conserving the total angular momentum. In fact, a star 
steadily becomes more centrally condensed and its radius is en- 
larged as it evolves from the ZAMS to the TAMS. This decreases 
the upper hmit of total angular momentum with time that the star 
can bear. 




J I I I I \ ^ I A I 



0.6 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
t/tMS 

Fig. 10. v/vzAMs ratio of equatorial rotational velocities as a function of 
the t/tus time ratio in the MS for stars with masses from M = 1.5 to 
3 Mq. vzams is the equatorial rotational velocity on the ZAMS, t is the 
stellar age and ?\is is the time that a star of given mass M can remain 
in the MS. Two extreme cases are represented: stars evolving as rigid 
rotators (instantaneous complete angular momentum redistribution) and 
as differential rotators without any angular momentum redistribution. 
Dashed lines represent the limits for J/M 0, while the solid lines are 
for stars that on the ZAMS have a specific angular momentum that in 
the TAMS becomes the critical one. 



5.1. Two limiting cases for ttie evolution of the equatorial 
rotational velocities 

Two simple and extreme cases of the evolution of the surface an- 
gular velocity can be considered, depending on the way the star 
redistributes its angular momentum during the MS evolutionary 
phase. For the limited values of rotational energy stored by the 
stars on the ZAMS stated above, / < /J^^, this can happen in 
two ways: (i) by entirely redistributing its angular momentum at 
any instant; (ii) by impeding any exchange of angular momen- 
tum among the stellar shells, so that each of them conserves its 
specific angular momentum. In the first case the star evolves as 
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Fig. 11. Theoretical evolution of the v/Vcnt velocity ratio calculated for 
several stellar masses, (a) Dependence with time of the v/vcrit ratio for 
several values of the total specific angular momentum 7, when stars 
evolve as rigid rotators all long the MS phase; (b) evolution of the ratio 
v/Vcrit normalized to its value on the ZAMS when rotators conserve the 
angular momentum by shells. 



a rigid rotator, while in the second case, we say the star evolves 
as a simple differential rotator. 



5.1 .1 . Stars evolving as rigid rotators 

The equatorial velocity at a time t elapsed from the ZAMS {t = 
0) in a star rotating as a rigid rotator is roughly given by 



v(0 
Vzams 



R(t) 



^ZAMS 



'ZAMS 



7(0 



(6) 



where [v(f), vzamsL WS), /?zams] and [7(f), 7zams] are the stellar 
equatorial linear rotational velocity, the radius and the moment 
of inertia at time I and at the ZAMS, respectively. We calcu- 
lated two-dimensional models of stars with rigid rotation as de- 
tailed in Appendix A. Figure 10 shows the evolution of the ratio 
v/vzAMS for stars with masses M - 1.5 and 3 Mq that conserve 
their total angular momentum during the MS evolutionary phase. 
The calculation was made for values of the angular momentum 
parametrized as / = ^ x Jj^^^, where q takes the values 0.1 
and 1. For objects undergoing total redistribution of the angular 
momentum, the velocity ratio v/vzams increases very slowly or 
remains nearly constant over the first two-thirds of the MS evo- 
lutionary period. Only in the last third of the MS phase do the 
changes of the stellar structure as a function of the angular mo- 
mentum produce some effects on the equatorial velocity, which 
slightly depend on the stellar mass. 
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The variation of the v/Vcdt ratio in rigid rotators of different 
masses and with total angular momenta J = qx J^^^ (q = 0.1, 
1/3, 2/3, 1) is plotted in Fig. 11a. The v/vcrit velocity ratio of rigid 
rotators remains nearly unchanged during the first two-thirds of 
the MS evolutionary phase, and increases sensitively only in the 
last third of the MS. However, (v/vcrit)zAMS clearly depends on 
the total angular momentum J. 



5.1 .2. Stars evolving as simple differential rotators 

The geometrical deformation of the stellar surface mainly de- 
pends on the specific angular momentum on the surface (Zorec 
et al., 2011a,b). Owing to the particular differential rotation as- 
sumed here, the specific angular momentum will become in- 
creasingly centrally condensed, so that the amount left in the 
surface will decrease strongly as the star evolves from ZAMS 
to TAMS. A simple account of the stellar deformation induced 
by the differential rotation could in principle be performed us- 
ing a homoeoidal description of the internal density distribution 
(Chandrasekhar, 1969), where the specific angular momentum 
is constant over spheroidal shells. However, owing to the fairly 
small effects carried by the rotation studied here on the stellar 
shape, we assumed that stars are spherical. The angular momen- 
tum of a spherical shell of width dr is 



dJ(t) 



f 



(7) 



where p(r) and Q(r) are the density and the angular velocity 
at radius r. 0(r) is assumed to be uniform over each shell, 
i.e. "shellular" distribution of the angular velocity. If r and r' 
are the radii at times t and t' of the same shell, whose mass 
dM(r) = 4npir)r^dr and angular momentum are conserved, we 
have 



Q'(r') 
0(0 



r) \r') ' 



(8) 



which for the equatorial velocity at t and f in the stellar surface 
implies that 



v{t) 



Reit'Y 



(9) 



where Rs(t) and Rs(t') are the equatorial radii of the star at times 
f and t', respectively. The change of the internal angular velocity 
as described by Eq. 8 induces a small variation on the moment of 
inertia that is taken into account as described in Appendix B . The 
variation of the equatorial velocity with time thus obtained and 
normalized to its value on the ZAMS is plotted in Fig. 10. The 
time dependent v/v'zams curves are also function of the stellar 
mass, while they are quite independent of the value of J. 

In Fig. 11 we show the variation of v/vcrit as a fimction of 
time for stars with masses from M = 1 .5 to 3 Mq, and several 
values of the angular momentum J - qx J^^^^ with q ranging 
from 0.1 to 1. Although the ratio v/vah depends on the value 
of J, its variation with time has a tiny dependence with mass. 
To demonstrate this small dependence with mass, the curves in 
Fig. lib are normalized to the respective values on the ZAMS: 
[v(M)/v(M)crit]zAMS- Owing to the low values of J chosen here, 
model stars can never have a critical specific angular momentum 
in the equator. For internal rotation laws given by Eq. 8 this can 
happen only for much higher values of J, which means that in 
this case the stars behave in the MS always as neat differential 
rotators. 




0.4 0.6 
t/tMS 

Fig. 12. Critical equatorial rotational velocities of stars with ditferent 
masses as a function of age in the MS. At each t/tus the objects are 
assumed to be rigid rotators. 



5.2. Observed evolution of rotational velocities vs. predicted 
extreme angular momentum redistribution regimes 

A general picture derived using observations of the evolution of 
rotational velocities in the MS as a function of mass is given 
in Fig. 13, where the average ratio of true rotational velocities 
(v'/vcrit) is shown as a function of the fractional age f/fMS in the 
MS for masses in the interval 1.7 < M < 3.3 Mg. Vcrit is the 
equatorial critical velocity given by 



: 436.7 



/ M V'^ 


\ 1 


1/2 


\Mol 


RcntiMj) 





kms \ 



(10) 



where R^niiM, f) is the stellar radius at critical rotation. For each 
star, the mass- and time-dependent critical radius Rai^^M, t) was 
determined using two-dimensional models of rigidly rotating 
stars, whose characteristics are given in Appendix A. Figure 12 
plots Vcrit(A/, t) against f/fMS for several stellar masses ranging 
fromM= 1.5 to 3 Mq. 

To obtain the diagram of average true rotational velocities 
(^/^crit) shown in Fig. 13, each observed vsin? parameter was 
divided by the corresponding v^aiM, t). The values were aver- 
aged in ruiming boxes with age bins of 0.18 for tjtMS < 0.3 and 
0. 15 for t/tus > 0.3 and mass bins of 0.5 Mq . These intervals are 
wide enough to ensure that the number of stars entering each av- 
erage leads to statistically reliable transformations from (v sin ;) 
to (v) and they are narrow enough to warrant that useful infor- 
mation on the evolutionary characteristics of rotation is hardly 
not obliterated. The true equatorial velocities come from the av- 
erage of the V sin / parameter corrected for sin / according to the 
statistical principles given by Chandrasekhar & Miinch (1950). 
The two-dimensional distribution was smoothed with a smooth- 
ing parameter derived as in Sect. 4 (see Bowman & Azzalini, 
1997). 

The two mass regions already identified in Fig. 7 depict in 
Fig. 13 an evolution of <v/Vcrit) with different characteristics. 
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Fig. 13. Distribution of the ratio v/vcit of true rotational velocities as a function of the fractional age in the MS evolutionary phase for normal stars 
in the mass interval 1.7 < M < 3.3 Mq. The filled circles are the average positions of the stars in running boxes with bins of 0.15 in f//MS and 0.5 
in M/Mq. 



The limit between these regions is roughly at M a; 2.5 Mq. To 
obtain a better insight on the characteristics of the evolution of 
the (v/vcrit) ratio in each of these mass-intervals, constant mass 
curves were extracted from Fig. 13 for three stellar masses: 2, 
2.5 and 3 Mq, which represent the average evolution of the ra- 
tio (v/v'ciit) as a function of time f/fMS and also describe the 
most typical aspects observed of the evolution of (v/vcrit)- These 
curves can be taken as tracers of the evolution of the average true 
rotational velocity against the fractional age f/fMS- However, for 
masses M > 2.5 Mq, these curves carry information on the evo- 
lution of two mixed groups of stars that we have akeady differ- 
entiated into slowly and rapidly rotating populations, of which 
the last is the most numerous one. Because the evolution of the 
rotational velocity in these stellar groups can present systematic 
differences, we aimed at obtaining the (v/vcrit) curves that rep- 
resent only the rapidly rotating stellar population. To this end, 
we used the correction factor (Table 6) that converts the mean 
values (v.s+r) of mixed slowly and rapidly rotating stars "s+r" to 
the searched rapid one (v,) = A:,(f) x (vs+i). The values calculated 
for the mid-time of each age interval were then interpolated for 
the required ages in the time interval < f/fMS ^ 1- Because 
Ki-(t) remains fairly constant over the first two thirds of the MS 
lifetime span, and the curves of average rotational velocities that 
we study here were normalized to the average velocity on the 
ZAMS, only the (v) velocities in the last evolutionary phases on 



the MS do undergo some little changes as compared to the orig- 
inal mixed average velocity distributions. 

The corrected (v/v'crit) curves as a function of the fractional 
age t/tus are presented in Fig. 14, where the vertical error bars 
give an insight on the uncertainties affecting the (v/vcrit) -ratio 
determination and the horizontal error bars are the averaged 
^t/iMs (from Table 1) in each bin. In this figure are superimposed 
the bands that coiTespond to the evolution of v/vait ratio for rigid 
rotators of masses from M = 1 .5 to 3 Mq previously shown 
in Fig. 11. They represent two angular momenta: J - J^^^^ 
(middle of the diagram) and J - (2/3)J^^^^ (lower right cor- 
ner). Even if in Fig. 14 the evolution of simple differential rota- 
tors is not shown, we notice that the evolution of (v/ v'ciit) differs 
strongly from what is predicted from the two extreme possibili- 
ties of angular momentum redistribution described in Sect. 5.1. 
In fact, in the first third of the MS and for both mass-groups, the 
ratio (v/vci-it) increases faster than suggested by the theoretical 
predictions. The (v/v'crit) ratio of stars with M - 2 Mq increases 
monotonically until the TAMS, while in stars with M > 2.5 Mq 
it reaches a maximum at tjtyis ~ 0.3. There is then a decrease 
that lasts roughly At/ t^s ~ 0.2, followed by a uniform value of 
<v/vcrit) until the TAMS. 

We calculated the mass- and time-dependent Va-it coiTe- 
sponding to all points in the curves of Fig. 14 and derive the 
evolution of the average true equatorial velocities (v) normalized 
to (v)zAMS for the same masses. The result is plotted in Fig. 15, 
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Fig. 14. Evolution of the equatorial velocity ratio (v/Vcnt) in the MS 
life time span. The ratios of true equatorial velocities are calculated for 
three masses: 2Mo, 2.5 Mq and SMq. The shaded regions correspond 
to the theoretical evolution of the v/Van ratio in rigid rotators plotted 
in Fig. 11, whose total angular momenta are / = Jj^^ (middle) and 
J = {2/3)J^^ (lower right comer). 




Fig. 15. Evolution of the equatorial velocity ratio (v/vzams) in the MS 
life time span, where Vzams is the true equatorial velocity of stars on the 
ZAMS. The average true equatorial velocities were calculated for the 
same masses as in Fig. 14. The shaded regions correspond to the theo- 
retical evolution of the v/vzams ratios shown in Fig. 10: rigid rotators 
(middle), simple differential rotators (lower left comer). 



where we also show the respective uncertainties. In this figure 
are superimposed the sequences of model v/v/ams variations for 
rigid (middle of the figure), and for simple differential rotators 
(lower left corner) previously shown in Fig. 10. 

The common characteristic of (v/vzams) curves seen in 
Fig. 15 is that for all masses they increase fast in the first third 
of the MS phase. For stars with M > 2.5 Mq (curves b and c), 
the ratio (v/vzams) decreases more or less monotonically in the 
second half of the MS phase. Owing to the shown uncertainties, 
we cannot be sure that the slopes reveal a slightly faster decrease 
than predicted for simple differential rotators. If it were actually 
the case, this decrease could imply some redistribution of angu- 
lar momentum toward the center of stars, provided that in the 
meantime no strong loss of angular momentum through stellar 
winds occurs. For stars with M = 2Mq (curve a), (v/vzams) re- 
mains almost constant from t/tus ~ 0-4 until t/tus ~ 0-7, and 
decreases slowly as the evolution approaches the end of the MS 
phase. This decrease follows the slope of simple differential rota- 
tors very closely, nearly as if the stars were not undergoing any 
redistribution of the angular momentum in the external layers 
during the final stages of their MS phase. 

In Figs. 14, 15, and 16 we show the uncertainties associated 
to the averages of fractional ages t/tus, which in all cases are 
'^i/'Ms ~ 0.1. In Table 2 we see that only for t/tus Z 0.5, sys- 
tematic deviations in the age estimations can be introduced by 
effects caused by the rapid rotation if Q/Qc ^ 0.8. However, 
because we note from Fig. 14 that (v/vcrit) ± o-v/vcm ^ 0.7 in 
t/tus ^ 0.5, they probably do not strongly affect the description 
of the evolution of rotational velocities obtained here. 

5.3. Observed evolution of rotational velocities vs. detailed 
calculations of the internal angular-momentum 
redistribution 

Ekstrom et al. (2008) computed stellar models to derive the evo- 
lution of surface velocities during the MS lifetime. Figure 16a 
displays the comparison of the observationally inferred evolu- 
tion of the v/Vcrit ratio with those obtained by Ekstrom et al. 
(2008) for three solar mass stars with different initial v/vcrit ra- 
tios. The theoretical v/v^it ratios reveal an initial short period 
where a star redistributes its internal angular momentum pass- 
ing from a rigidly rotating object to a differential rotator at the 
low-energy regime (Maeder & Meynet, 2000). We must note 
that the observed curve corresponds to a mixed stellar popula- 
tion where the average mass is M = 3 Mq. The averaging of 
velocities is necessarily made over a distribution of initial ve- 
locities that we do not know. We could assume some kind of 
MaxwelUan or Gaussian initial distribution of equatorial veloc- 
ities to obtain curves for single initial v/Vcrit ratios using per- 
haps some kind of deconvolution. The operation would never- 
theless lead to uncertain results. However, the differences shown 
in Fig. 16a are eloquent enough by themselves to indicate that in 
the first half of the MS evolutionary phase the angular momen- 
tum in actual stars may undergo other redistribution processes 
than those theoretically predicted today. This difference appears 
stiU more clearly in Fig. I6b, where we compare the observed 
evolution of v/vzams ratios with theoretical curves for the same 
initial velocity ratios as in Fig. 16a. In this figure the incidence 
of initial velocity distributions would have hardly any effect in 
the first half of the MS evolutionary phase, so that the above 
noted discrepancy between theory and observation seems to be 
confirmed. 

Excepting the very first drop of velocity ratios to be caused 
by an initial fast internal angular redistribution, we note that in 
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Fig. 16. Comparison of the observed evolution of equatorial velocity 
ratios in the MS life time span for stars with M = 3 Mq inferred from 
observations with theoretical ones calculated by Ekstrom et al. (2008) 
for different v/vcnt on the ZAMS. (a) Evolution of (v/vciit) ratios; (b) 
evolution of <v/vzams) ratios. 



the time interval 0.1 < t/tus ^ 0.6 the theoretical curves in 
Fig. 16b differ very Uttle from the evolution of rigid rotators. In 
the last MS evolutionary span, the detailed calculations suggest a 
lowering of the equatorial velocity whose steepness is flatter than 
for simple differential rotators, which indicates that the stellar 
surface layer gains some angular momentum from that stored in 
the stellar interior. 

In all the above comparisons we note that for stars with 
masses M > 2.5 Mq there seems to exist a characteristic step 
or time scale at which significant changes of (v) are produced 
that do not seem to be present in the theoretical models, which 
amounts to 



6t ~ 0.2 tMS- 



(11) 



5.4. Discussion of tine possible internal distribution of the 
anguiar veiocity 

A different insight in the evolution of rotational velocities can 
be obtained by comparing the minimum total specific angular 
momentum J/M required to account for the equatorial veloci- 
ties shown in Fig. 15. The minimum angular momentum J/M is 
meant here to be the value stored in a rigid rotator that accounts 
for the observed rotational velocity v at a given t/tMS- Let us then 
assume that the sequences of (v/vzams) plotted in Fig. 15 repre- 
sent the evolution of the rotational velocity in actual single stars 
with the masses indicated in the figure. Because the mass loss 
rate in these stars is very low, we may consider that they evolve 
with conservation of the total angular momentum all along their 
MS phase. For all masses in Fig. 15 and for the following stages 
f/fMS ~ 0.07 (near the ZAMS), t/tus ~ 0.3-0.4 (when v is max- 
imum for stars with M > 2.5 Mq), t/tus ~ 0.65 (near maxima 
and/or inflections of v) and t/tms ~ 0.9 (near the TAMS), we 
calculated the total specific angular momentum J/M that rigid 
rotators require to account for the observed equatorial velocities. 
We also estimated the angular momentum of the critical rigid ro- 
tators at each chosen evolutionary stage. The estimated values of 
J/M are given in Table 7 and the results show that 

(i) for all evolutionary stages and in all studied masses we 
have on average J/M < (J/M)^^^, but roughly J/M > 
0.5 X (7/M)ZAMS). 

(ii) for almost all masses and evolutionary stages ranging from 
t/tMS = 0.3 to 0.6, we notice that J/M > (//M)zams. 
which indicates that to account for the rotational veloc- 
ity of stars at these stages, more rotational kinetic energy 
is needed than we inferred when we assume that they are 
mere rigid rotators on the ZAMS; 

(iii) for masses M = 2 Mq we find that J/M < (J/M)1^^^ in 
all evolutionary stages, while for higher masses stars need 
J/M > (J/M)J^^ to account for their surface rotation at 
t/tus ~ 0.3. 

We can then comment that 

(iv) the condition J/M < (J/M)^^^^ in all evolutionary stages 
means that the evolution of rotational velocities in the 

MS of stars in the studied mass ranges can be consistent 
with low regimes of rotational energies, i.e. t = K/\W\ < 
r(ZAMS)ff; 

(v) knowing that for all evolutionary stages (including the 
ZAMS) the total specific angular momentum is estimated 
supposing that stars are rigid rotators, values J{t)/M > 
(y/M)zAMS at any t/tyis > imply that there must be 
some transfer of angular momentum from the center to- 
ward the surface. However, because rigid rotation in the 
initial evolutionary phases demands that Qcore = ^^enveiope> 
in latter evolutionary phases where J/M > (J/M)zams 
would necessarily imply Qcore < ^^enveiope. which could be 
not realistic. An increase of i^enveiope could then be possi- 
ble if stars actually evolve as differential rotators having 
Qcore > ^^enveiope since the ZAMS. However, the condition 
J/M < (J/M)^^^ suggests that they could be simple dif- 
ferential rotators on the ZAMS; 

(vi) because in the studied range of stellar masses the total 
angular momentum of stars is conserved during the MS 
evolutionary phase, in the stages where it is inferred that 
J/M > (J/M)^^^^, stars have to redistribute their total an- 
gular momentum and end their MS phase behaving as neat 
differential rotators. 
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Table 7. Synopsis of the average evolution of equatorial velocities of 
2 Mo, 2.5 Mo and 3 M^ stars. 
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Noting from the above arguments that stars can behave as 
differential rotators, we can try to speculate 

(A) what the internal rotation law of a star on the ZAMS can 
look Uke, if we assume that its total specific angular mo- 
mentum is the highest we can find in a given MS evolution- 
ary stage as shown in Table 7; 

(B) what the required values of JjM and the associated inter- 
nal rotation laws can be to account for the observed aver- 
age rotational velocity at the spotted evolutionary stages in 
Table 7, if the star is assumed to be a differential rotator, 
and what the possible internal rotation laws on the ZAMS 
are that can afford the same JjM values to explain the aver- 
age ZAMS equatorial velocities. 

The ratio f2o/^e strongly depends on the characteristics of the 
internal rotational law (Qq is the angular velocity of the center 
of the star, r - RIRe - 0; Q.^ is the angular velocity in the 
equator at the surface, r - RIR^ = 1). These notations should 
not be confused with Qcore and i^enveiope. which vaguely refer to 
angular velocities in larger domains in the core and in the stellar 
envelope, respectively. Because at the moment we do not know 
what this law can be, we adopted the model-laws inferred by 
Maeder & Meynet (2000) for early-type stars in the MS to make 
an educated guess on the possible ratios ilo/Q.^. These laws can 
be sketched analytically as follows (Zorec et al., 2007) 

Q(r) = Qo[l -px exp(-a/r^)], (12) 

where Q(r) is uniform over each spherical shell of radius r, 
i.e. "shellular" internal distribution of the angular velocity; from 
Eq. 12 it is = ^(r = !)• The quantity p is a contrast parame- 
ter; a is a parameter that depends on the radius of the stellar core, 
which is assumed to be rotating rigidly; b describes the steepness 
of the change of il(r) from Q.o to £2e. We note that p - Ois suited 
for rigid rotation and p - 1 leads to the strongest possible Qo/f2e 
ratio for Rayleigh's stability condition d jidr > to be satisfied. 
The values of b range from b « 3.5 on the ZAMS to > 2 
in the TAMS (Maeder & Meynet, 2000). To make quantitative 
estimates, let us adopt a 2.5 Mg test star with the rotational char- 
acteristics displayed in Table 7. The results obtained in the frame 
of the above quoted hypotheses (A) and (B) are given in Table 8. 

According to question (A), we see that for the 2.5 Mq star 
the highest JjM value in the MS occurs at f/fMS = 0.36, where a 
rigid rotator needs iJIM)t/i^^=036 = 1-27 x 10^^ cm^ s"', which 



is also the minimum value possible for (JIM) to account for 
the observed average equatorial velocity v = 231kms~^ at 
this evolutionary stage. Assuming that we have (//M)zams = 
(7/M),/,„5=()36, using the rotation law in Eq. 12 we obtain the 
contrast parameter pzAMS - 0.25 and see that (Qo/i^e)zAMS = 
1.31 to explain v = 182 kms"' on the ZAMS. 

Following the question/hypothesis (B), we ascribed several 
values to p and asked which 7/M values explain v = 231 km s"' 
at f/fMS = 0.3, if the 2.5 M© test star rotates with the law given 
by Eq. 12. The obtained values 7/M = f(p) and ratios (Qo/^^e) 
are given in the upper part of Table 8 tmder the label "(B)". We 
then required that the />zams values according to the 7/M = f{p) 
values obtained can account for the specific rotational velocity 
V = 182 km s"' on the ZAMS (lower part of Table 8 under the la- 
bel "(B)"). From Table 8 it is clear that the higher the value of the 
contrast parameter p, the higher JjM needs to be to explain the 
same average equatorial velocity. We also notice that as soon as 
/?ZAMS > 0.53, //M exceeds (JIMf^^^^ = 1.97 x lO'^cm^s^i, 
so that the star must be a neat differential rotator already on 
the ZAMS. In Fig. 17 we show the rotation laws sketched with 
Eq. 12 that would correspond to the parameters p and (Qo/^^e) 
given in Table 8 for the ZAMS and for tjtus = 0.36. 



Table 8. Contrast parameters p, angular velocity ratios flo/i^e and total 
specific angular momenta JIM in a 2.5 Mo test star. 



f/fMS = 0.36 
P 

(7/M)/10"cm2s-' 


Questions/hypotheses 
(A) (B) 


0.00 0.20 0.40 0.60 0.80 1.00 
1.00 1.21 1.52 2.06 3.19 7.02 
1.27 1.48 1.81 2.36 3.52 7.47 


f/fMS = 0.00 

PZAMS 

(flo/i^e)zAMS 

(7/M)/10"cm2s-' 


(A) (B) 


0.25 0.39 0.53 0.67 0.81 0.96 
1.31 1.56 1.96 2.63 4.03 8.80 
1.27 1.48 1.81 2.36 3.52 7.47 



6. Comments and conclusions 

Paper III provides rotational velocity distributions for A-type 
star groups defined according to their spectral types. The deriva- 
tion of fundamental parameters allowed us to disentangle the ef- 
fects of mass and evolutionary stage in these distributions. In 
spite of this more refined classification, the striking features al- 
ready observed in Paper III remain unchanged (Figs. 6 and 9): 
(i) the noticeable lack of slow rotators in the low-mass end, i.e. 
M < 2.5Mq (ii) the genuine bimodal distributions for the high- 
mass end (M > 2.5Mq). Taking advantage of the results shown 
in Figs. 13, 14, and 15, we can attempt to put forward the fol- 
lowing comments on some properties of these distributions. 

On the one hand, we saw in Figs. 14 and 15 that in interme- 
diate and late A-type stars (1 .6 < M/Mq < 2.5) the rotational ve- 
locities undergo an acceleration from tltus ~ to titus ~ 0.3- 
0.4 and then remain high for a fairly long time, which can cer- 
tainly help to keep the fraction of slow rotators low during sub- 
sequent evolutionary stages up to the TAMS. 

On the other hand. Fig. 15 shows that stars with masses 
2.5 < M/Mq < 3.5 also accelerate their surface rotation in the 
time interval ranging from f/fMS ~ to f/fjvis ~ 0.3, but undergo 
an efficient deceleration up to the TAMS, which we already com- 
mented in Sects. 5.2 and 5.3. We could then think of an initial. 
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Fig. 17. Internal rotation law.s Q.(r)/n„ sketched with Eq. 12 and in- 
ferred for a 2.5 Mq test star at t/tus = 0.36 and on the ZAMS, according 
to several parametrized values of the contrast parameter p adopted for 
t/tus - 0.36 that all account for the observed rotational velocity at this 
evolutionary phase. For each parameter Pt/»Ms we obtained the corre- 
sponding pzAMs, so that for the same specific total angular momentum 
J/M the corresponding angular velocity ratio flo/^\ can explain the 
observed rotational velocity vzams given in Table 8. 



more or less genuine, bimodal distribution of rotational veloci- 
ties among stars with masses in the 2.5 < M/Mq < 3 interval. 
The commented deceleration of rotational velocities can help to 
renew the population of slowly rotating stars as the evolution 
proceeds. 

We conclude, however, that we have no substantial evidence 
to say why stars in the low-mass interval are more or less system- 
atically not subject to efficient enough decelerating mechanisms 
that enable them to start their MS evolutionary phase with fairly 
high rotational velocities. In the same way the fraction of ZAMS 
slow rotators remains unexplained as well. We can guess, nev- 
ertheless, at two possibilities: (i) some stars in the present high- 
mass interval could start their MS phase after having completely 
exhausted their angular momentum in the pre-MS phase through 
magnetic braking, since it is difficult to believe that they never 
had any angular momentum (Tassoul, 1978); (ii) there could be 
a fraction of undetected binaries where the tidal braking could 
have some effect. A systematic search of the fraction of rapidly 
rotating stars among those with small vsinj parameters could 
perhaps shed new light on this question. 

We have studied the evolution of the ratio of average true ro- 
tational velocities (v/vcnt) (v is the true equatorial velocity; Vcnt 
is the critical equatorial velocity) of MS stars with masses rang- 
ing from 1.7 Mq to 3.3 Mq. Using 2D models of rotating stars we 
calculated the Vcrit for individual masses and ages, and obtained 
the evolution of the average true equatorial velocity (v) for three 
stellar masses: 2, 2.5 and 3 Mq plotted in Figs. 14 and 15. These 
masses summarize the main characteristics of the evolution of 
rotational velocities in the MS of stars in the 1.7 < M/M© < 3.3 
mass range. 



The observed evolution of rotational velocities strongly dif- 
fers from that theoretically predicted for two limiting cases of 
internal angular momentum redistribution: (i) rigid rotation (to- 
tal redistribution); (ii) conservation of the specific angular mo- 
mentum by each stellar shell (no redistribution). For all studied 
masses, the observed trends suggest that some effective mecha- 
nism of angular momentum redistribution favors the acceleration 
of the rotation in the stellar surface during the first third of the 
MS phase. This also makes stars in the studied mass-range pos- 
sibly behave as differential rotators during all their MS phase. 
After the acceleration in the first third of the MS phase, stars 
with masses M = 2 Mq evolve with little change of their equa- 
torial rotational velocity in the second third of the MS. In the 
last third of the MS, the rotational velocities decelerate because 
they were conserving the specific angular momentum per shell 
without any redistribution. Stars with masses M > 2.5 Mq seem 
to have an efficient deceleration of the surface rotation during 
the second third of the MS, which is faster than imposed by a 
regime of no redistribution of angular momentum. A new ac- 
celeration of the surface rotational velocity seems to happen in 
the 0.6 < t/tus ^ 0.7 time interval, but it decelerates again in 
the last third of the MS phase. This might be caused by some 
redistribution of angular momentum toward the center of stars, 
since according to their low mass-loss rates they cannot undergo 
strong losses of angular momentum. The variation of the surface 
rotational velocities proceeds at characteristic time scales that 
are mass-dependent, i.e. 6t « 0.2 ?ms- 

The comparison of the observed evolution of rotational ve- 
locities for stars with M = 3 Mq with that predicted by detailed 
calculations of redistribution processes in the stellar interiors 
shown in Fig. 16 reveals very strong deviations between theory 
and observations in the first half of the MS evolutionary phase. 
We note that a more pronounced spin-down trend during the sec- 
ond part of the MS compared with the M = 3 Mq model from 
Ekstrom et al. (2008) is found by HGM. According to the present 
work, this faster spin-down seems to be present in the last third 
of the MS, as shown in Fig 16b. 

From the total specific angular momentum calculated to ac- 
count for the observed velocities of stars in different evolution- 
ary stages if ffiey were rigid rotators, we conclude that stars of all 
masses studied here could start the MS evolutionary phase with 
a lower total angular momentum than the limiting value ascribed 
to critical rigid rotators on the ZAMS. Also, the observed behav- 
ior of ffie equatorial velocities suggests that they evolve during 
most of the MS phase as simple differential rotators. However, 
the angular momenta obtained for masses M > 2.5 Mq in their 
last third of the MS phase suggest that they must have larger 
amounts of rotational energy than rigid critical rotators can bear 
in the TAMS. This indicates that by the end of the MS they 
evolve as neat differential rotators. 

Using a test star with M = 2.5 Mq, we have speculated on 
the properties that its internal rotation may have on the ZAMS, 
if at this stage the surface rotational velocity corresponded to the 
total angular momentum required at t/tus - 0.36, which is the 
highest estimated for stars of this mass and in particular, it is 
higher than the expected one on the ZAMS behaving as rigid ro- 
tator. We conclude that an angular velocity ratio ilo/Q.^ up to 1.3 
on the ZAMS could be expected. Assuming that the 2.5 Mq test 
star is a simple differential rotator (J/M < (J/M)'^^^), its equato- 
rial velocity can correspond to angular velocity ratios Oo/Qg as 
high as 2.6 on the ZAMS, so that Qo/Oe ^ 2.1 for p = 0.6 
at t/tus = 0.36. If the object were a neat differential rotator 
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{JIM > {JIMfJ^), stable OJQe ratios up to 8.8 could exist 
on the ZAMS and up to 7 at tltus = 0.36. 

The present discussion assumes that the internal rotation of 
stars has a shellular-like distribution. We noted in Paper 111, how- 
ever, that some A-type stars could present a radiative/convective 
structural dichotomy that could induce other types of angular 
velocity laws in the external layers and could accordingly af- 
fect the interpretation of the measured rotational velocities. This 
dichotomy might be responsible for anomalous von Zeipel co- 
efficients in fast rotating A-type stars (Zhao et al., 2009). Fast 
rotation can, in particular, affect the temperature gradient in the 
stellar envelopes, producing thus enlarged convective zones in a 
similar way as was suggested for massive fast rotators (Maeder 
et al., 2008). Convection and rotation may then be coupled some- 
how to induce other rotational laws than of mere shellular type 
(Zorec et al., 2011a,b). In the future, tests should then be made 
to see which internal rotation profiles can actually prevail. These 
tests could perhaps be carried out by studying the non-radial pul- 
sation modes of early A-type stars. 

Acknowledgements. We are thankful to an anonymous referee for her/his valu- 
able comments and suggestions that helped to improve the presentation of our 
results. We are also very grateful to Prof Hugo Levato for the communication of 
his data. We thank Mrs A. Peter for her prompt and efficient language editing of 
this paper. This work has widely made use of the CDS database. 



References 

Abt, H. A. 2009, AJ, 138, 28 

Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359, (ALG) 
Abt, H. A. & Morrell, N. 1. 1995, ApJS, 99, 135 
Attridge, J. M. & Herbst, W. 1992, ApJ Lett., 398, L61 
Barnes, S. A. 2003, ApJ Lett., 586, L145 

Bodenheimer, P. 1981, in Fundamental Problems in the Theory of Stellar 
Evolution, lAU Symp. 93 (D. Reidel Publishing Co., Dordrecht), 5 

Bowman, A. W. & Azzalini, A. 1997, Applied smoothing techniques for data 
analysis: the kernel approach with S-plus illustrations, Oxford statistical sci- 
ence series No. 1 8 (Clarendon Press, Oxford) 

Castelli, F. 1991, A&A, 251, 106 

Castelli, F. & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, lAU 
Symp. 210, ed. N. Piskunov, W. W. Weiss, & D. F. Gray (ASP, San Francisco), 

20P 

Chandrasekhar, S. 1969, Ellipsoidal figures of equilibrium, ed. Chandrasekhar, 

S. (Yale University Press, New Haven) 
Chandrasekhar, S. & Miinch, G. 1950, ApJ, 111, 142 
Choi, P 1. & Herbst, W. 1996, AJ, 1 1 1, 283 
Clement, M. J. 1974, ApJ, 194, 709 

Debemardi, Y. 2000, in Birth and Evolution of Binary Stars, lAU Symp. 200, ed. 

B. Reipurth & H. Zinnecker (ASP, San Francisco), 161 
Deutsch, A. J. 1970, in Stellar rotation, lAU Colloquium No.4, ed. A. Slettebak 

(D. Reidel Publishing Co., Dordrecht), 207 
Diaz, C. G., Gonzalez, J. E, Levato, H., & Grosso, M. 201 1, A&A, 531, A143 
Edwards, S., Strom, S. E., Hartigan, P, et al. 1993, AJ, 106, 372 
Ekstrom, S., Georgy, C, Eggenberger, P., et al. 2011, ArXiv e-prints 
Ekstrom, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 
Endal, A. S. & Sofia, S. 1981, ApJ, 243, 625 
Erspamer, D. & North, P 2003, A&A, 398, 1 121 
ESA. 1997, The Hipparcos and Tycho Catalogues, ESA-SP 1200 
Fekel, F C. 2003, PASP 1 15, 807 

Fremat, Y. & Zorec, J. 2003, in Modelling of Stellar Atmospheres, lAU Symp. 

210, ed. N. Piskunov, W. W. Weiss, & D. F Gray (ASP San Francisco), 27P 
Fremat, Y, Zorec, J., Hubert, A.-M., & Floquet, M. 2005, A&A, 440, 305 
Fremat, Y, Zorec, J., Hubert, A. M., Floquet, M., & Cidale, L. 2003, in 

Modelling of Stellar Atmospheres, lAU Symp. 210, ed. N. Piskunov, 

W. W. Weiss, & D. F Gray (ASP San Francisco), 26P 
Guthrie, B. N. G. 1982, MNRAS, 198, 795 
Hauck, B. & Mermilliod, M. 1998, A&AS, 129, 431 

Herbst, W., Bailer-Jones, C. A. L., & Mundt, R. 2001, ApJ Lett., 554, L197 
Howarth, I. D. 2004, in Stellar Rotation, lAU Symp. 215, ed. A. Maeder & 

P. Eenens (ASP, San Francisco), 33 
Huang, W., Gies, D. R., & McSwain, M. V. 2010, ApJ, 722, 605, (HGM) 
Hubrig, S., North, P, & Mathys, G. 2000, ApJ, 539, 352 



Jefferys, W. H., Fitzpatrick, M. J., & McArthur, B. E. 1998a, Celest. Mech., 41, 

39 

Jefierys, W. H., Fitzpatrick, M. J., McArthur, B. E., & McCartney, J. E. 1998b, 
GaussFit: A System for least squares and robust estimation. User's Manual, 
Dept. of Astronomy and McDonald Observatory, Austin, Texas 

Keppens, R., MacGregor, K. B., & Charbonneau, P 1995, A&A, 294, 469 

Konigl, A. 1991, ApJ Lett., 370, L39 

Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Introduction to Stellar Winds, 
ed. Lamers, H. J. G. L. M. & Cassinelli, J. P. (Cambridge University Press) 

Landstreet, J. D., Bagnulo, S., Andretta, V., et al. 2007, A&A, 470, 685 

Lang, K. R. 1992, Astrophysical Data I. Planets and Stars., ed. Lang, K. R. 
(Springer- Verlag, Berlin) 

Lanz, T. & Catala, C. 1992, A&A, 257, 663 

Latour, J., Toomre, J., & Zahn, J.-P 1981, ApJ, 248, 1081 

Levato, H. & Grosso, M. 2004, in Stellar Rotation, lAU Symp. 215, ed. 
A. Maeder & P. Eenens (ASP, San Francisco), 51, (EG) 

Lucy, L. B. 1974, AJ, 79, 745 

MacGregor, K. B. & Brenner, M. 1991, ApJ, 376, 204 

Maeder, A., Georgy, C, & Meynet, G. 2008, A&A, 479, L37 

Maeder, A. & Meynet, G. 2000, ARA&A, 38, 143 

Matsumoto, T, Hanawa, T, & Nakamura, F. 1997, ApJ, 478, 569 

Moon, T. T. & Dworetsky, M. M. 1985, MNRAS, 217, 305 

Mouschovias, T. C. & Morton, S. A. 1985a, ApJ, 298, 205 

Mouschovias, T. C. & Morton, S. A. 1985b, ApJ, 298, 190 

Napiwotzki, R., Schoenberner, D., & Wenske, V. 1993, A&A, 268, 653 

Pinsonneault, M. H., Kawaler, S. D., & Demarque, P 1990, ApJS, 74, 501 

Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, A&A, 424, 727 

Pudritz, R. E. 1985, ApJ, 293, 216 

Reiners, A. & Royer, F 2004, A&A, 415, 325 

Renson, P & Manfroid, J. 2009, A&A, 498, 961 

Richardson, W. H. 1972, J. Opt. Soc. Am., 62, 55 

Royer, F. 2009, in Lecture Notes in Physics, Vol. 765, The Rotation of Sun and 
Stars, ed. J.-P. Rozelot & C. Neiner (Springer Verlag, Berlin), 207 

Royer, F, Gerbaldi, M., Faraggiana, R., & Gomez, A. E. 2002a, A&A, 381, 105, 
(Paper I) 

Royer, F, Grenier, S., Baylac, M.-O., Gomez, A. E., & Zorec, J. 2002b, A&A, 

393, 897, (Paper 11) 
Royer, F, Zorec, J., & Gomez, A. E. 2007, A&A, 463, 671, (Paper 111) 
Sackmann, 1. J. 1970, A&A, 8, 76 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 
Schatzman, E. 1962, Ann. Ap., 25, 18 

Sheather, S. J. & Jones, M. C. 1991, J. R. Statist. Soc B., 53, 683 

Shu, R H., Najita, J. R., Shang, H., & Li, Z.-Y. 2000, in Protostars and Planets 

IV, ed. V. Mannings, A. P. Boss, & S. S. Russell (University of Arizona Press, 

Tucson), 789 

Silverman, B. W. 1986, Density estimation for statistics and data analysis. 
Monographs on Statistics and Applied ProbabiUty No. 26 (Chapman & Hall, 

London) 

Slettebak, A., Colhns, I. G. W, Boyce, R B., White, N. M., & Parkinson, T. D. 

1975, ApJS, 29, 137 
Soderblom, D. R., Jones, B. F, & Fischer, D. 2001, ApJ, 563, 334 
Soderblom, D. R., Stauffer, J. R., MacGregor, K. B., & Jones, B. F. 1993, ApJ, 

409, 624 

Spruit, H. C. 1999, A&A, 349, 189 
Spruit, H. C. 2002, A&A, 381, 923 

Stauffer, J. R., Hartmann, L., Soderblom, D. R., & Bumham, N. 1984, ApJ, 280, 
202 

Stjpieh, K. 2000, A&A, 353, 227 

Tassoul, J.-L. 1978, Theory of rotating stars (Princeton University Press) 

van Leeuwen, F. 2007, A&A, 474, 653 
von Zeipel, H. 1924, MNRAS, 84, 665 

Zhao, M., Monnier, J. D., Pedretti, E., et al. 2009, ApJ, 701, 209 
Zorec, J. & Briot, D. 1997, A&A, 318, 443 
Zorec, J., Cidale, L., Arias, M. L., et al. 2009, A&A, 501, 297 
Zorec, J., Fremat, Y, & Cidale, L. 2005, A&A, 441, 235 

Zorec, J., Fremat, Y, Delaa, O., et al. 201 la, in Active OB stars: structure, evo- 
lution, mass loss, and critical limits, lAU Symp. 272, ed. C. Neiner, G. Wade, 
G. Meynet, & G. Peters (Cambridge University Press), 103-105 

Zorec, J., Fremat, Y, & Domiciano de Souza, A. 2007, in ASP Conf. Ser., Vol. 
361, Active OB-Stars: Laboratories for Stellar and Circumstellar Physics, ed. 
A. T. Okazaki, S. P Owocki, & S. Stefl (ASP San Francisco), 542 

Zorec, J., Fremat, Y, Domiciano de Souza, A., et al. 201 lb, A&A, 526, A87 

Zorec, J., Fremat, Y, Hubert, A. M., Floquet, M., & Cidale, L. 2002, in Rev. Mex. 
Astron. Astrofis. Conf. Ser., ed. J. J. Claria, D. Garcia Lambas, & H. Levato, 
Vol. 14, 114 



18 



Zorec & Royer: Rotational velocities of A-type stars. IV., Online Material p 1 



Appendix A: Models of rigid rotators 

Only the general dynamical aspects induced by rotation were 
considered here to calculate the mass distribution in a star and 
consequently the gravitational potential of a centrifugally dis- 
torted star. We separated the primary dynamical effects produced 
by rotation from those induced by evolution. The primary ther- 
modynamic effects carried by the stellar evolution were taken 
into account using barotropic relations calculated with stellar 
models without rotation. We assumed therefore that the changes 
produced in the P - P(p) relation at a given evolutionary stage 
of a star by the several instabilities and the diffusion of chemi- 
cal elements unleashed through the stellar evolution by rotation 
have second-order effects on the establishment of the dynamical 
equilibrium of the rotating star. In principle, one could use the 
barotropic relations derived with models of stellar evolution with 
shellular rotation, but the results will not be more rehable. This 
approach is used and discussed recently by Zorec et al. (2011b) 
for massive and intermediate-mass fast rotating stars. 

Because we are not interested in the precise description 
of all non-linear time-dependent phenomena associated with 
the viscosity and with internal flows in rotating stars, and be- 
cause the total energy carried by the meridional circulation is 
low, our models are axisymmetric, steady state and circula- 
tion free. Because of these assumptions, our model-stars behave 
as barotropes [Poincare-Wavre theorem (Tassoul, 1978)]. We 
adopted internal rotational laws of conservative form, O = 
where m is the distance to the rotation axis. The rigid rotation 
is a special case of this type of rotational laws. In this case, 
the gravitational potential <J)(zzt,z) and the density distribution 
pinr, z) in the rotating star are simultaneous solutions to the hy- 
drostatic equihbrium equation: 

p-iyP = W^ + fvj-^e^, (A.l) 
and to the Poisson equation 

AO) = AnGp, (A.2) 

where (m, (p, z) are the cylindrical coordinates with z containing 
the rotation axis; is the unit vector perpendicular to the z-axis; 
P is the pressure; j = D.w^ is the specific angular momentum. 

Equations (A.l) and (A.2) are solved with the adopted com- 
plementary barotropic relation 

P(p) = ap^" + bp^" , (A.3) 

where the constants a, b, ja and jb were adjusted to (i) repro- 
duce the pressure Pc and the density pc in the center of the non- 
rotating star of given mass and evolutionary stage; (ii) ensure a 
continuous distribution of the pressure-density relation at the ra- 
dius of the stellar core; (iii) obtain the right stellar mass at the 
stellar radius as tabulated by Schaller et al. (1992) for ID evolu- 
tionary models for the initial metalUcity Z = 0.02. The function 
(A.3) is continued in the stellar atmosphere by another pressure- 
density relation calculated by Castelli & Kurucz (2003) for stel- 
lar atmospheres as a function of the parameters (Teff , log g). 

The first-order effects from the stellar evolution are thus ac- 
counted for by the pressure-density relations in the center of the 
star and by the dP/dp gradients. An additional term in relation 
(A.3) could in principle also take into account the presence of 
the convective regions in the envelope induced by fast rotation, 
but we did not do this here. The only rotational effect considered 
here on the P = P(p) relation is through the mass-compensation 
effect (Sackmann, 1970), which increases the density p^ in the 
center of the star. For this, we iterated pc until the nominal stellar 



mass M was obtained. This iteration also impUes that the central 
pressure Pc changes in accordance. 

The gravitational potential 3>(ziT, z) was obtained by solving 
Poisson equation (A.2) with the cell-method adapted by Clement 
(1974) for stellar structure calculations. The density distribution 
p(cT, z) was derived through the integrated form of (A.l). 

Given an angular rotational velocity Q. and a barotropic re- 
lation (A.3), the solutions of equations (A.l) and (A.2) were 
performed over the entire space. The iteration of O and p was 
stopped when the highest density difference in the {m, z)-space 
is max(5p/p) 10"^. In our iterations the virial relation 6 = 
[2{K+U)-W]/\W\ = OiK = kinetic energy; U = internal energy; 
W = total gravitational potential energy) is verified to better than 
(5 » 2 X 10-"* in the ZAMS models and (5 =^ 6 x 10"^ by the TAMS 
models. Since in the frame of conservative rotational laws the 
surfaces of constant pressure, density, and of total potential are 
parallel, the rotationally distorted shape of our models is defined 
by the total equipotential surface that contains the polar "photo- 
spheric" radius Rp. This radius is identified by the layer whose 
density satisfies the model-atmosphere relation tross(p) = 2/3 in 
the stellar atmosphere models of Castelli & Kurucz (2003). The 
local effective temperature at the pole needs to also satisfy the 
gravity darkening effect. We accordingly modified the effective 
temperature given by Schaller et al. (1992) for the given mass M 
using von Zeipel's approximation (von Zeipel, 1924). The trans- 
formation to the rotation dependent effective temperature was 
performed following the procedure given in Fremat et al. (2005). 

The models calculated in this work are for specific total 
angular momentum scaled as J/M = 0.1, 1/3, 2/3 and X 
(y/M)J^^, where iJ/M)^^^^ is the critical value supported by 
the stars on the TAMS. We also assumed that the stars evolve 
conserving their initial angular momentum. The velocity ratios 
for J/M = and 1 X (y/M)TAMS plotted in Fig. 10. The 
characteristic parameters of the models calculated are given in 
Table B.l. 

Appendix B: Models of differential rotators 

The angular velocity distribution Q'(r') at time f can be derived 
at time t if we obtain the relation r' - r'(r) that describes the 
change of the radius r of a shell at time t to the radius r' of the 
same shell at time t'. This can be done observing that the mass 
inside a given radius r at time t must be the same as that inside 
r' at time f: 

f pix')x'^dx = f pix)x^dx. (B.l) 
Jo Jo 

A similar relation can be obtained from (7) by imposing that 
the angular momentum inside r at time t is the same as inside 
r' at the time t' . This relation is much more difficult to man- 
age, however, since the distribution of the angular velocity £l(x') 
and f must be iterated. In Fig. B.l we show the internal rota- 
tional profiles at different evolutionary phases identified by the 
ratios f/fMS^ where t is the stellar age from the ZAMS and tus 
is the stellar age on the TAMS. These distributions were ob- 
tained assuming that on the ZAMS the star is a rigid rotator. 
In Fig. B.l the angular velocities are normaUzed to the angular 
velocity at r = 0. The shaded region corresponds to the extent of 
the stellar core. We note that in current models of massive and 
intermediate-mass rotating stars, it is assumed that in the con- 
vective core Ocore - constant. In that case a distribution like the 
one sketched for t/tus = 0.247 is expected. 

To include the effect of the internal rotation on the stellar 
moment of inertia and on the stellar radius, we used the angular 
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velocities shown in Fig. B.l that are interpreted as conservative, 
i.e. written as a function of the distance ■ui to the rotational axis. 
Since the total rotational energy is low, i.e. JIM<{JIM)cnu dif- 
ferences in the estimation of stellar radii are small as noted by 
the loci of points in Fig. 10 for J/M = and J/M-{J/M)cnt. 




Fig. B.l. Internal angular velocity at different evolutionary stages in 
MS for a 2 Mq rigidly rotating star on the ZAMS that evolves without 
any exchange of the angular momentum among the stellar shells. The 
shaded zone represents the stellar core, where the current models with 
rotation assume that it rotates rigidly. In this case, the internal angular 
velocity is redistributed as sketched for the t/tus -0.247 phase. 
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Table B.l. Parameters characterizing the models of rigid rotators that evolve with constant angular momentum J=qx J^^, with for rigid 
rotation. 
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Table B.l. continued. 
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Notes. For each mass and for ^ = 0.1 is given the Vcnt, while for other p values is given the ratio v/Vcnt- f/fms is the fractional steUar age, fjvis being 
the time the star of mass M can spend in the MS. Q. is the angular velocity and Slcnt is the critical angular velocity. Pc is the density at the stellar 
center RpIRo is the polar radius in solar units and R^IRp is the equatorial to polar radii ratio, v is the equatorial rotational velocity and Vcrit the 
equatorial critical rotational velocity. KI\W\ is the ratio of kinetic to gravitational potential energy. 



